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The low-energy QCD predictions to be tested by the DIRAC experiment are 
revised. The experimental method, the setup characteristics and capabilities, 
along with first experimental results are reported. Preliminary analysis shows 
good detector performance: alignment error via A mass measurement m.\ = 
1115.6 MeV/c^ with cr = 0.92 MeV/c^, pjT relative momentum resolution 
CTQ f» 2.7 MeV/c, and evidence for vr+vr^ low momentum Coulomb correlation. 
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1 Introduction 



Quantum Chromodynamics (QCD), responsible for the strong interaction sec- 
tor of the Standard Model (SM) has successfully been tested only in the per- 
turbative region of high momentum transfer (Q > 1 GeV) or at short relative 
distance Ar ^ Ti/Q (Ar < 0.2/m). Here the constituent quarks behave as 
weakly interacting, nearly massless particles. The QCD in the perturbative 
region, as any gauge theory with massless fermions, presents chiral symmetry. 

In the nonperturbative region of low momentum transfer (low-energy), 
say Q < 100 MeV, or equivalently at large distance (Ar > 2/m), asymptotic 
freedom is absent, and quark confinement takes place. In the low energy 
region the chiral symmetry of QCD must be suontaneously broken. 

The Chiral Perturbation Theory (ChPT)Er0 seems to be the candidate 
theory for low energy processes. It exploits the mechanism of spontaneous 
breakdown of chiral symmetry (SBChS), or, in other words, the existence of a 
quark condensate. In order to test the existence of the quark condensate, the 
particularly significant symmetry breaking effect refers to the 5*- wave tttt scat- 
tering lengths. From the theoretical point of view, tttt scatteringJa^ a deeply 
studied problem. Within standard ChPT, Gasser and LeutwylertloEl and also 
Bijnens and coUaboratorsD as well as within ther,GjeneraJized Chiral Perturba- 
tion Theory (GChPT), Stern and collaboratorsH'B, have obtained expressions 
for the TTTT scattering amplitude in the chiral expansion. j.. 

The leading order expansion of the scattering amplitude ig3 

Ais; t,u)=a^ + ^i^s- ^M|) + (1) 

where a and /3 encode the information on the strength of the quark condensate. 
In the limit of a strong quark condensate, one has a « 1, /? w 1. A substantial 
departure of a from unity signals a much smaller value of the condensate. 

The values predicted for the isospin / = and 1 = 2 5-wave scattering 
lengths flg and a§ can be confronted with the future expierimental values of 
the DIRAC experiment. The available experimental dataQ for the scattering 
length is Og = 0.26 ± 0.05. Based on these data there is no possibility to 
estimate the strength of the quark condensate and so to measure the extent 
of chiral symmetry breaking— 

The DIRAC experimenlH aims to determine the difference of the scat- 
tering lengths A = Ioq — OqI with 5% accuracy, by measuring the lifetime 
of pionium (tt+tt^ bound state). For the first time experimental evidence in 
favour of or against the existence of a strong quark condensate in the QCD 
vacuum could be within reach. 
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2 Pionium lifetime 

Pionium is a metastable bound state, produced by tt"*" and 7r~ electromagnetic 
interaction and decaying into tt'^tt'^ due to strong interaction. To obtain pion 
scattering lengths in a model independent way, a measuremeot of the hfetime 
of pionium has been proposed many years ago by NemenovB. The measure- 
ment of the pionium lifetime (r) will allow to determine the difference Ioq — csqI 
of the strong 5- wave TTTr-scattering lengths for isospin 7 = and 1 = 2. 
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Figure 1. The pionium lifetime, in units of 10^ s, as a function of the combination 
(ajj — ciq) of the 5- wave scattering lengths. The band delineated by the dotted lines takes 
account of the uncertainties , coming from th eoretical evaluations, low energy constants and 



The general expression for the pionium decay width (F), according to 
ChPThEl at leading and next-to-leading order in isospin breaking, is 

l=T=l-a'-p*-{4-al + e)'il + K), (2) 

where a is the fine structure constant ; 

p* = (m2^ - aPo - a^Ml^/Af'^ is the CM momentum of vr" ; 
e = (0.58'± O.ief - 10-2 ; \ ^ 1.07 • lO'^. 
Using Eq.(|), = 0.206 and af, = -0.0443, Gasser et al.tahave evaluated 
the pionium lifetime in the ground state r = 3.25 • 10~^^s. In the jsospin 
symmetry limit e = 0; A' = 0, Eq.(^) becomes the Deser type formulauJ. 

Eq.(||) allows to get a relative error of the scattering lengths difference 

^/^f'"!} « 5%, if DIRAC measures the hfetime with a ^ w 10% error. 
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According to GChPTy the T vs. {qq — Gq) dependence is presented in 
Fig.|l|. Hence the experimental determination of the pionium hfetime could 
be interpreted in quark condensation terms for three different cases: First, 
if the central value of r is close to 3 x 10"^^ s, lying above 2.9 x 10^^^ s, 
then the strong condensate assumption of ChPT is firmly confirmed, since its 
predictions of (oq — Uq) lie between 0.250 and 0.258. Second, if the central 
value of T lies below 2.4 x 10^^^ s, it is the scheme of GChPT, which is 
confirmed, since the corresponding central value of a would lie above 2. The 
third possibility is the most difficult to interpret. If the central value of r 
lies in the interval 2.4 4- 2.9 x 10~^^ s, then, because of the uncertainties, 
ambiguities in the interpretation may arise. 



3 Experimental method 

In an abundent production of oppositely charged pions the Coulomb interac- 
tion can form atomic tt+tt^ bound states. If the pions have a small relative 
momentum in their CM system (Q ~ 1 MeV/c) and are much closer than the 
Bohr radius (387 fm), then pionium atoms are produced with a high produc- 
tion probability due to the large wave function overlap. Such pions originate 
from short-lived sources (p, ct;,A), and not from long-lived ones (rj^Kg), be- 
cause in the latter case the separation of the two pions is in most cases larger 
than the Bohr radius. 

The pionium production cross section is proportional to the double in- 
clusive cross section da^ / {dpidp2) for tt'^tt" piairs from short-lived sources, 
without Coulomb interaction in the final statecl, and to the squared atomic 
wave function of nS'-states at the origin |'(/'„(0)|^ 

^ = (2.)3|i|^„(0)r-^, ^ (3) 

dp A Ma dpidp2 U=p^=EA 

where pA , Ea and Ma are momentum, energy and mass of the pionium atom 
in the Lab system respectively; pi and p2 are the 7r+ and tt^ momenta in the 
Lab system, and they must satisfy the relation pi — P2 — Pa/'2 to form the 
atomic bound state. Atoms formed in this way are in a S'-state. 

After production in hadron-nuclcus interaction relativistic pionium atoms 
(2 GeV/c < PA < G GeV/c) are moving in the target. They can decay or, due 
to electromagnetic interaction with the target material, get excited or broken- 
up (ionized). Using atomic interaction cross sections, for a given target ma- 
terial and thickness, one can calculate the break-up probability for arbitrary 
values of pionium momentum and lifetime. In Fig.H there are presented these 
dependencies for the pionium momentum p = 4.7 GeV/c. 
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Figure 2. Probability of pionium break-up in the target. 

Comparison of the measured break-up probability P^r = ua/Na (ratio 
of broken-up - and produced - pionium atoms) with the calculated 
dependence of Pbr on r (see Fig.^ gives a value of the lifetime. 

The break-up process gives characteristic tt+tt^ pairs, called atomic pairs. 
They have a small relative momentum in their CM system (Q < 3 MeV/c), a 
small opening angle {6± « 6/7 « 0.35 mrad for pA = 4.7 GeV/c) and nearly 
identical energies in the Lab system {E+ = E- at 0.3 % level). 

3.1 Determination of broken-up pionium atoms (ua) 

The measurement of broken-up ua pionium atoms is realized through the 
analysis of the experimental distribution in Q of tt+tt" pairs. 

The free pion pair distribution can be written as the sum of the non- 
Coulomb (nC) (no final state interaction) and the Coulomb (C) pair distribu- 
tion (Fig|): 




Figure 3. Accidental, non-Coulomb (long-lived sources) and Coulomb (short-lived sources) 
pion pair production and pionium production 
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dQ ~ dQ ^ dQ 
dm^ dN^^P dN" 



acc 



; - HQ)Ac{Q){l + aQ) 



dQ dQ ' dQ 

where Ac{Q) is the Coulomb and (1 + aQ) is the strong correlation factors. 

Here we assumed that the non-Coulomb distribution of tt+tt" pairs (with- 
out FSI) can be extracted from the experimental distribution of accidental 
pairs '^{Q). The free pion pair distribution is then given by 

+ ^om) [f + MQKl + aQ)] , (4) 

A'o, /, a - free parameters. In the region Q > 3 MeV/c, there are mainly free 
pairs. This part of the distribution is fitted with the function (^) containing 
the experimental distribution of tt+tt" accidental pairs. The extrapo- 

lation of the approximation function to the region Q < 2 MeV/c yields the 
number of free pairs in this region. Hence the value ua of atomic pairs is 

Q<2 

From the measured break-up probability Pbr = nA/NA^ where Na is the 
calculated total number (according to (|^)) of produced pionium atoms, and 
from the dependence of Pbr on the lifetime r, one can derive a pionium ground 
state lifetime and hence a value for A = Ioq — AqI. 



4 Experimental setup 

The experimental setuplll (Fig.^ has been designed to detect pion pairs and 
to select atomic pairs at low relative momentum with a resolution better than 
1 MeV/c. It was installed and commissioned in 1998 at the ZT8 beam area 
of the PS East Hall at CERN. After a calibration run in 1998, DIRAC has 
been collecting data since summer 1999. 

The 24 GeV/c proton beam extracted from PS is focused on the target. 
The secondary particle channel, with an aperture of 1.2 msr, has the reaction 
plane tilted upwards at 5.7° relative to the horizontal plane. It consists of 
the following components: 4 planes of Micro Strip Gas Chambers (MSGC) 
with 4 X 512 channels; 2 planes Scintillation Fiber Detector (SciFi) with 2 x 
240 channels; 2 planes Ionization Hodoscope (IH) with 2 x 16 channels; 1 
Spectrometer Magnet of 2.3 Tm bending power. Downstream to the magnet 
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Figure 4. Experimental setup 

the setup splits into two arms placed at ±19°, relative to the central axis. 
Each arm is equipped with a set of identical detectors: 4 Drift Chambers 
(DC), the first one common to both arms and with 6 planes and 800 channels, 
the other DCs have altogether per arm 8 planes and 608 channels; 1 Vertical 
scintillation Hodoscope (VH) plane with 18 channels; 1 Horizontal scintillation 
Hodoscope (HH) plane with 16 channels; 1 Chcrcnkov detectors (Ch) with 10 
channels; 1 Preshower scintillation detector (PSh) plane with 8 channels; 1 
Muon counter (Mu) plane with (28+8) channels. 

For suppressing the large background rate a multik'vcl trigger was de- 
signed to select atomic pion pairs. The trigger levels are defined as fol- 
lows: To = (VH ■ PSh)i ■ {VH ■ PSK)2 ■ IH, fast zero level trigger; Ti = 
{VH ■ HH ■ Ch ■ PSh)i ■ {VH ■ HH ■ Ch ■ PSh)2, first level trigger from the 
downstream detectors; T2 = Tq ■ {IH ■ SciFi), second level trigger from the 
upstream detectors, which selects particle pairs with small relative distance; 
Tz is a logical trigger which applies a cut to the relative momentum of parti- 
cle pairs. It handles the patterns of VH and IH detectors. T3 did not so far 
trigger the DAQ system, but its decisions were recorded. 

An incoming flux of ^ 10"'^^ protons/s would produce a rate of secondaries 
of about 3 X 10^/s in the upstream detectors and 1.5 x 10^/s in the downstream 
detectors. At the trigger level this rate is reduced to about 2 x lO'^/s, with an 
average event size of about 0.75 Kbytes. 

With the 95 /xm thin Ni target, the expected average pionium yield within 
the setup acceptance is ~ 0.7 x 10~^/s, equivalent to a total number of ~ 10^^ 
protons on target to produce one pionium atom. 

5 First experimental results 

The data taking has been done mainly with tt+tt" and pTr~ pairs and also e+e" 
pairs for detector calibration. For the first data analysis only the most simple 
events were selected and processed, those with a single track in each arm. 
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with signals in DC, VH and HH. The tracks in the DCs were extrapolated to 
the target plane crossing point of the proton beam. A cut was applied along 
X and Y distances between the extrapolated track and the hit fiber of the 
SciFi planes (ISmm divided by the particle momentum in GeV/c, to take into 
account the multiple scattering effect). Finally, these events were interpreted 
as TT+TT" or pi:~ pairs produced in the target. 

The difference in the time-of-flight Ai between the positive particle (left 
arm) and negative particle (right arm) of the pair at the level of VH is pre- 
sented in Fig]|. 




Figure 5. VH time-of-fiight difference distribution Figure 6. Positive particle momen- 
for pair events turn versus VH time-of-flight dif- 

ference for particle pairs 

The first interval —20 < < —0.5 ns corresponds to accidental hadron 
pairs (mainly tt+tt^). In the second interval —0.5 < At < 0.5 ns one observes 
the peak of coincidence hits associated to correlated hadron pairs over the 
background of accidental pairs. The width of the correlated pair peak yields 
the time resolution of the VH {at ~ 250 ps). The asymmetry on the right side 
of the peak is due to admixture of protons in the 7r+ sample, that are p7r~ 
events. Hence the third interval 0.5 < At < 20 ns contains both accidental 
pairs and p-K^ events. 

This time-of-flight discrimination between 7r"'"7r~ and p7r~ events is effec- 
tive for momenta of positive particles below 4.5 GeV/c. This is demonstrated 
in Fig.^, where the scatter plot of positive particle momentum versus dif- 
ference in time-of-flight At in VH is shown. The single particle momentum 
interval accepted by spectrometer is 1.3 7.0 GeV/c. 

For correlated tt+tt" pairs Coulomb interaction in the final state has to 
be considerated, because it increases noticeably the yield of tt+tt^ pairs with 
low relative momentum in CM (Q < 5 MeV/c). For accidental pairs this 
enhancement is absent. 
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Figure 7. Correlated it~^tt pairs with po- Figure 8. pn invariant mass for proton mo- 
sitive particle momenta piab < 4.5 GeV/c menta piab > 3 GeV/c. 
and Qt < 4 MeV/c. 



Fig.0 shows the distribution of the longitudinal component Ql (the pro- 
jection of Q along the total momentum of the pair) for correlated pairs. There 
are plotted pair events with positive particle momentum piat < 4.5 GeV/c, 
occuring within the "correlated" At peak and with transversal component 
Qt < 4 MeV/ c to increase the fraction of low relative momentum pairs. 

In the region \Ql\ < 10 MeV/c there is a noticeable enhancement of 
correlated tt+tt^ pairs due to Coulomb attraction in the final state. 

The most important parameter for data analysis is the resolution in 
and Qt- This has been measured by the reconstruction of the invariant mass 
of pn^ pairs. The distribution of pvr^ invariant mass is presented in Fig.^. 
Positive particles are restricted to momenta larger than 3 GeV/c, and the 
time-of-flight must lie in 0.5 < At < 18 ns. A clear peak at the A mass 
TOA = 1115.6 MeV/c^ with a standard deviation a ~ 0.92 MeV/c^ can be 
seen. These mass parameter values show a good detector calibration and 
coordinate detector alignment, with an accuracy in momentum reconstruction 
better than 0.5 % in the kinematic range of A decay products. This gives for 
the relative momentum resolution trg ~ 2.7 MeV/c. For tt+tt^ pairs a better 
resolution can be obtained, due to the different kinematics. 

6 Conclusion 



The DIRAC setup test and calibration have been done, and the DIRAC expe- 
riment began data taking. To achieve the goal - measurement of the pionium 
lifetime with 10% precision - we have to consider a number of at least 20000 
recorded " atomic pairs" . Improvements in hardware and software will con- 
tinue this year. These will result in a better data quality. 
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